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Abstract 

In Hofava-Lifshitz gravity, the extra dynamical scalar mode could play significant role in cos- 
mology. However it has been pointed out that such a scalar may suffer from the strong coupling 
problem in IR. We address this issue in this paper. Our analysis shows that the scalar mode could 
decouple naturally at A = 1 due to the extra gauge symmetry. On the other hand, the fact that 
the scalar mode becomes ghost when 1/3 < A < 1 is a real challenge to the theory. We try to 
overcome this problem by modifying the action such that the RG flow lies outside the problematic 
region. We discuss the cosmological implications of the action and calculate the power spectra of 
scalar and tensor modes. 
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I. INTRODUCTION 



Diffeomorphism is essential to Einstein's relativity theory of gravity It has been widely 
believed to be exact in any theory of gravity. However, in the recent proposal by Hof ava[j], y] 
on gravity theory, it is no longer an exact symmetry. The basic idea behind Hofava's theory 
is that time and space may have different dynamical scaling in UV limit. This was inspired 
by the development in quantum critical phenomena in condensed matter physics, with the 
typical model being Lifshitz scalar field theory 0, 0]. In this Hofava-Lifshitz theory, the time 
and space will take different scaling behavior as 

x -> 6x, t -> b z t, (I) 

where z is the dynamical critical exponent characterizing the anisotropy between space and 
time. Due to the anisotropy, instead of diffeomorphism, we have the so-called foliation- 
preserving diffeomorphism. The transformation is now just 

t -> t(t) 

x i -> ?(x j ,t). (2) 

As the result of this "reduced" gauge symmetry, there are more physical degrees of 
freedom in the theory. In fact, there exist an extra dynamical scalar degree of freedom in 
Hofava-Lifshitz gravity, as shown in [H, 0]. This scalar deg ree of freedom and its physical 
implication has also been discussed in cosmology @, 3l[. However, it was pointed out 



in 0] that this scalar mode would be strongly coupled to the matter in the IR fixed point 
A = 1, at which one would expect the recovery of diffeomorphism. If this is the case, it 
would lead to unacceptable effects in many experiments. In the next section, we would like 
to show that this would not happen. 

The key point in our analysis is that we take the point of view that the diffeomorphism is 
only an approximate symmetry even at IR. In fact, as we will review shortly, it is not hard 
to see that even at IR fixed point, there exist various other terms, involving spatial higher- 
derivative terms, which break the diffeomorphism, even though they should be very much 
suppressed. This is very different from the case in Fierz-Pauli's massive gravity In the 
massive gravity theory, there exist extra physical degree of freedom. It was a serious issue on 
how this degree of freedom get decoupled in the massless limit, where the diffeomorphism is 
completely recovered^, lo| . In our case, we will show manifestly that the extra scalar degree 



of freedom could be decoupled without trouble, due to the existence of extra gauge symmetry 
at IR fixed point rather than the complete recovery of diffeomorphism. The breakdown of 
full diffeomorphism at IR fixed point also suggest that the usual Stuckelberg trick could not 
be used directly, especially taking into account of the projectability condition. 

Another issue on the scalar mode is whether it is a real physical degree of freedom. 
The debate in the literature focus on if one should choose the lapse function to be only 
the function of time, or in other words, if the lapse function should be projectable. The 
different choice seems lead to completely different physics. For example, it was found that 
without the projectability condition there were new static spherically symmetric solutions 
to Hofava-Lifshitz gravity and its modifications [111 1~2|. These new solutions may have 
profound physical implications in solar system tests HJ • However, it was proved in fl3[ that 
these new solutions do not respect the projectability condition. From our point of view, 
taking the lapse function as the function of time is the most natural choice. With this 
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choice, the gauge transformation looks transparent and simpler. Moreover the Hamiltonian 
constraints form a closed algebra, and the theory gets rid of the pathology found in [l|, [TEJ . 
Simply speaking, the theory is well-defined with the projectability condition. As a result, 
the extra scalar mode becomes the physical one and the key ingredient in our following 
discussion. 

On the other hand, this extra scalar mode could not be always physical. When the pa- 
rameter A lies between 1/3 and 1, this mode is actually a ghost. We would like to emphasize 
that the existence of the ghost is the real challenge to the original Hofava-Lifshitz gravity 
theory. It indicates that the theory may neither well-defined at UV, nor UV complete. In 
particular, at UV, one wish that the theory becomes non-relativistic and the speed of light is 
much larger than the constant one at IR. In the original proposal of Hofava-Lifshitz gravity, 
this requires that the theory stay near A = 1/3. However in IR, one may expect that the 
theory would flow to A = 1. As the RG flow is from A = l/3toA = l, the theory inevitably 
suffers from the existence of ghost. To get away from this trouble, we try to modify the 
action in a way that the RG flow may be from UV with A > 1 to IR with A = 1. To simplify 
the analysis, instead of considering the most general form of the action, we only consider 
the potential with the marginal terms and the most relevant terms. This will be the topic 
in section 3. 

In the remaining part of this paper, we discuss the cosmological implications of the 
scalar mode, and also calculate the power spectra of scalar and tensor perturbations. The 



cosmology of the Hofava-Lifshitz gravity has first been discussed in [111. 1 171 Il8l. Il9l|. and then 
widely studied in the literature 20] from various angles. In this paper, we take the extra 
scalar as an alternative to the inflaton and study its power spectrum. We also calculate 
the power spectrum of tensor mode in modified Hofava-Lifshitz gravity action proposed in 
section 3. In our treatment, we simply ignore the RG flow and use the standard technology 
in cosmological perturbation theory. The problem turns out to be quite similar to the trans- 
Planckian problem in inflation. Instead of the WKB approximation used in 0, 3(|, we 
apply the technology in trans-Planckian physics and study the equation of motion of scalar 
perturbation stage by stage. We find that the power spectra are scalar invariant. This is not 
a surprise since the classical evolution is a pure de-Sitter phase, which has time translation 
invariance. We also notice that the tensor-to-scalar ratio is sensitive to the time of horizon- 
crossing of tensor and scalar modes, and can be small if at the time of scalar crossing the 
horizon A is near 1. 

The paper is organized as follows. In section 2, we study the gravitational scalar in the 
Hofava-Lifshitz gravity theory. In section 3, we present our modification of the Hofava- 
Lifshitz gravity action. In section 4 and 5, we calculate the power spectra of the scalar and 
tensor perturbations respectively. We end with some discussions in section 6. 



II. SCALAR MODE IN HORAVA-LIFSHITZ GRAVITY 

Since time direction plays a privileged role in the whole construction, it is more convenient 
to work with ADM metric 

ds 2 = -N 2 dt 2 + g ij (dx i + N*dt) (dx j + N j dt). (3) 

Due to the anisotropy between time and space, the usual diffeomorphisms reduce to the 
foliation-preserving diffeomorphisms, generated by infinitesimal transformation: 

6t = S°(t), 5x l = C(t,x). (4) 
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The essential point is that £° is just the function of t. This leads to the following transfor- 
mations on the metric components: 



Sgij = di£ k g jk + dj£ k g ik + £ k d k gij + £°gij, 

= d i i j N j + edjNi + e 9lJ + i°Ni + e ^, 

SN = £ j djN + i°N + £°A>. (5) 

The above transformations could be obtained by taking a nonrelativistic limit of usual 
relativistic diffeomorphisms. It is more convenient and natural to choose N being just the 
function of t. There are a few advantages to work with this choice. With this choice, the 
gauge symmetry is simpler and transparent. Furthermore, in the Hamiltonian formulation, 
the constraints could form a closed algebra since the momentum conjugate to N does not 
lead to a local constraint [l( . As a result of less constraints than standard GR, the physical 
degrees of freedom in the theory include not only the massless gravitons but also another 
propagating scalar. The existence of extra scalar field has profound meaning in cosmology. 
In [6J, we showed that for the action without the detailed balance condition, this scalar may 
lead to scale invariant spectrum. 

On the other hand, if one abandon the projectability condition and let N be the function 
of both t and x\ one will find that the theory would be ill-defined, as shown in [l], [l5| . 

At the special value A = 1, the theory develops an enhanced time-independent U(l) gauge 
symmetry acting via 

6N i = die, 6 9ij =0. (6) 

Due to the existence of extra gauge symmetry, the scalar mode is not physical anymore. It is 
remarkable that even with this extra gauge symmetry, the total gauge symmetries is different 
from the usual diffeomorphisms in general relativity. In other words, the diffeomorphisms 
has not been recovered at A = 1. This fact is essential to understand why at A = 1 the extra 
scalar degree of freedom could be decoupled without trouble. 

To understand the decoupling better, let us consider the perturbation around FRW met- 
ric: 

ds 2 = —dt 2 + a{t) 2 5ijdx l dx^ 

= a 2 (r])(-dr] 2 + 5 ij dx i dx j ). (7) 

Here for simplicity we focus on the flat universe, and use the co-moving time rj = J dt/a 
as time variable. The above metric could be reduced to Minkovski spacetime if a(t) is a 
constant. Perturb the flat metric, and use the ADM formulism in co-moving time, 

ds 2 = -(M 2 - MU^drf + 2^4^ + g^dx^ (8) 

The fluctuations around the above metric could be 



N = a(n)(l + A), (9) 
Mi = afrXdiB + Vi), (10) 

9a = « 2 (^){(i - 2V0*j - d A E - 2d a F j) + hij}, (ii) 

where A,B,ip,E are scalar perturbations, Vi and Fj are vector perturbations, and is 
the gauge-invariant tensor perturbation describing gravitational wave. Under the gauge 
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transformations (jSJ), we have 



A 

B 
E 

V5 



A = A--{?a)', 

B = B — a(', 
E = E + 2(, 

a 



V; 



Vi + (&j 
hi 



where we have decompose the spatial vector as 

e = ei + 9x 



(12) 

(13) 
(14) 

(15) 

(16) 
(17) 
(18) 

(19) 



with being divergenceless and £ being a scalar. The gauge invariant variables besides hij 
are 



* = S + f, 
■V, = Vi - R. 



(20) 

(21) 
(22) 



We will only focus on the scalar perturbations. It is convenient to work with the gauge 

A = 0, E = 0. (23) 

Note that the above gauge choice is consistent with the projectability condition. Since £° 
is the only function of t, the gauge transformation on A would not spoil the projectability 
condition. 

If the scale factor is a constant, the above gauge transformations reduce to the ones in 
flat spacetime. In [l6j], the gauge invariant perturbations about the flat spacetime have been 
analyzed carefully. Actually, from the discussion there, one can see that the extra dynamical 
scalar mode can decouple without trouble. It was claimed in 1161 that such scalar mode is 
not a propagating mode. This is not true. The problem comes from the fact that the lapse 
function is projectable so that it induce a non-local super-Hamiltonian constraint. For the 
flat spacetime, the perturbation A of the lapse function is a pure gauge and can be set to 
zero safely. Even if A is kept nonvanishing, the variation with respect to A would only lead 
to non-local constraint, which is less powerful than the local one. 

In terms of ADM metric, the action of original Hofava-Lifshitz gravity theory can be 
written as[2| 



,s'„ / dtd 3 Xy/gN \ -^KijG ijkl Ku - 



K 

~2 



— -Ca — I R; 



-Rgij + A w g,ij 



Qijkl 



(24) 
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where Ky is the extrinsic curvature of the spatial hypersurface; Cy is the Cotton tensor which 
can be used to preserve the detailed-balanced condition in constructing the action; G ljkl is 
the De Witt metric on the space of metrics that preserve the anisotropic diffeomorphism, 
and Rij is the Ricci tensor in spatial hypersurface. Their definitions are 

Kij = ^(9ij - ViJV,- - VjNi), (25) 
Ca = e ikl V k (r^ - ±R5j) , (26) 
G ijU = 1 y.k g3 i + g ii g kij _ Xgijg u (27) 

Here and throughout the paper, a dot over the quantity means taking the derivative with 
respect to cosmic time t, while a prime denotes that to co-moving time rj. The first term in 
( 124)) involving only the extrinsic curvature is the kinetic term, while the others are potential 
terms. A is the coupling constant in the kinetic term, and runs expectedly to A = 1 at IR 
regime at which the kinetic term goes back to the one in the general relativity. This specific 
form of the action is governed by the detailed-balance condition, which is just applied by 
Hofavafor convenience to decease the number of arbitrary parameters. The expansion of the 
action gives 

f 2 K 2 K 2 II 

dtd 3 x^N I -ziKuir* - \K 2 ) + —h-e^RuVjrtk 

vy \k 2K 1 ' 2oo 4 3 2u 2 3 

~^f^R ij + (^R 2 + ^R - 3A^) } • (28) 

Comparing this action with the Einstein-Hilbert action in IR limit 

Seh = j d 4 x^N{(K tJ K^ -K 2 )+R- 2A}, (29) 

with x° = ct, we can recover the speed of light, Newton constant and the cosmological 
constant by the parameters introduced before, 



c=^/^, G=J*-, A = -A W . (30) 
4 V 1 - 3A' 32ttc' 2 w K J 

Thus at IR the theory recovers nearly the usual general relativity, with the higher derivative 
terms of spatial metric components as the modifications. Even though the higher derivative 
terms are highly suppressed at IR, strictly speaking, the theory always breaks diffeomor- 
phism, or locally Lorentz invariance. 

In Hofava's original paper ^ the coupling constant A runs to 1 in IR limit. And in UV, 
because of the anisotropy between space and time, the speed of light is not a constant and 



may be extremely large, which could be used to explain the horizon and flatness problem [18 
But from f)30p we know that this can only occur in the case A < 1/3 if we take A to be positive, 
taking into account of the fact A is directly related to cosmological constant. However, this 
raise the worry that the marginal coupling constant A can never run to its infrared value 
A = 1, which is directly in contrast with our former description. To solve this problem, it 



was proposed that one should do analytical continuation on the parameters 11 



fj, — > ifM, uj — > — iu, (31) 
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which leaves the action real. And under this continuation, we see from (1301) that 



A 



4 V3A-1' (32) 

and there is no conflict between A > and A > 1/3. And when A — > 1/3 proposed by 
Hofava as the ultraviolet value of this coupling constant, we have a very large speed of light, 
which can naturally solve the casuality problem in cosmology without inflation. 

To understand the extra scalar degree of freedom, let us come back to the scalar pertur- 
bation we studied before [1, 0]. We need only focus on the kinetic term in the above action, 
which is just 



S K = / dtd 3 x{ 3aa 3 (l - 3A) 



%b 2 

+ 6Hi[)ip + 9H 2 ip 2 



31-A 



(33) 



Several remarks are in order: 



1. From the action, it is obvious that the scalar mode ip is physical when A < 1/3 and 
A > 1, while when 1/3 < A < 1 the mode is a ghost, indicating that the theory is not 



well-defined [2 II]. At the special value A = 1, the mode is decoupled, as we will clarify 
more below. And at A = 1/3, the theory has extra symmetry, as discussed carefully 
in j2| • This fact is the same as the one found in [l], [2| where the perturbations around 
the flat spacetime were studied. 



2. More interestingly, the equation of motion of ip takes the following form: 

1 - 3A 



1- A 



+ .... 



(34) 



This indicates that when A — > 1, the scalar field ip could be decoupled naturally, in 
contract with the claim in [3]. It seems that the strong coupling problem does not 
exist in our case. 

3. The absence of the strong coupling problem may stem from the fact that we take 
different points of view on gauge transformations. In our case, we stick to the require- 
ment that the lapse function should be projectable, as originally advocated in [2|. As 
a result, we do not expect that the diffeomorphism is recovered at A = 1. Instead, 
the decoupling of the extra scalar mode comes from the fact that there is extra gauge 
symmetry at A = 1. This is conceptually different from the case studied in [7J and 
Fierz-Pauli massive gravity [8|. 

4. Technically it is remarkable the equation of motion of ip has a pref actor proportional 
to 1/(1 — A) rather than (1 — A). This difference has significant physical implication. 
In our case, this means that the scalar mode could be decoupled without trouble. 
Another way to see this is to cast the scalar mode into canonical form such that the 
mode become non-physical at A = 1. It is remarkable that in [H, 0], the equation of 
motion of the scalar mode around the flat spacetime background has the prefactor 
(1 — A). However this is due to different gauge choice. It has been shown in [l6[ by 
rescaling the field, one has the same equation of motion. In fact, no matter what kind 
of gauge choice, the physical dispersion relation is exactly the same. This suggests 
that for the cosmological perturbations, the different gauge choice would not lead to 
different dispersion relation. Namely, the extra scalar mode may decouple naturally 
as A — > 1. 
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III. MODIFIED HORAVA-LIFSHITZ GRAVITY 



The existence of the ghost is fatal to the theory. It means that the theory is not well- 
defined, not mentioning UV completeness. One may expect that we can always work in the 
region outside A G [1/3,1]. However this cannot be guaranteed, considering our ignorance 
of the details of RG flow. On the other hand, in the practical application in cosmology, one 
wish the RG flow is from A~l/3toA=lin original Hofava-Lifshitz gravity. In this 
paper, we take a modest attitude and try to modify the Hofava-Lifshitz gravity such that 
the RG flow may happen always with A > 1. In order to do so, we have to abandon the 
detailed balance condition. As it is well-known, the detailed balance condition may not be 



essential to the theory [26l. |28| . [33J. The imposing of such condition is pragmatic to simplify 
the action. In principle, one may relax this condition and consider more general form of 
the action. In this paper, we do not want to consider the most general form of the action. 
Instead, we just consider the marginal spatial kinetic part and most relevant deformations, 
besides the time kinetic terms. The action we start with is of the form 

S g = J dtd 3 x^N {a(K i:j K ij - XK 2 ) + £(X)R + <x(A) 

-/3 (/3 1 V i R jk V i H ih + foV % R jk Vm ik + foViKV'R) } . (35) 

Here we only keep the marginal terms that are power-counting renormalizable and domi- 
nant in UV limit, besides the lower-dimensional terms to recover IR behaviors. 1 The other 
marginal terms being cubic of Ricci scalar and Ricci tensor, and the other relevant terms like 
R 2 and RVR are neglected for simplicity. For a complete discussion on all possible terms 
maintaining the power-counting renormalizability, see [26| . 

Because of the breakdown of the detailed balance condition, the coupling constants before 
each terms are independent. The couplings could be connected to the speed of light, the 
Newtonian coupling constant and the cosmological constant of Einstein's general relativity 
in IR limit, 

c 2 = i (36) 
a 

16ttG = — , (37) 

ca 

A = (38) 

Here we see that c 2 can be positive, if we choose a proper form of the function £(A). Fur- 
thermore, we can require c to be very large when A is near its ultraviolet value. In Hofava's 
original paper, he suggested A — > 1/3 at the UV limit, which gives a large speed of light 
in (I30I) or (132]) . Here we only take this condition as a constraint on the function £(A). For 
instance if the theory requires A to be larger than the unity at UV as we will propose as a 
condition to exclude the ghost field, the function £(A) may be divergent when A tends to be 
infinity. 



1 Actually, after some integrals by parts and using the Bianchi identity, the /32V 'i-RjfcV J ' R tk term can be 
converted to a (/?2/4) Vji?V l i? term and some higher order terms. In our current work we can just set 
$3 — 03 + @2 /4 and discard the 02 term. This will change in the calculation of non-Gaussianity. Thanks 
Shinji Mukohyama for useful discussions. 
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For our use let us have a glance of the classical dynamics of the universe under such an 
action. In a homogenous and isotropic universe, 



ds 2 



-di 2 



a h'ij dx doc ■ 



h 



'.J 



Sa + 



L-/Cx 2 



(39) 



where /C is the parameter to describe the spatial curvature. Under this metric, the universe is 
homogeneous and isotropic, which will greatly simplify our following calculations. To apply 
our foliated diffeomorphism, we need use the ADM formalism of this Robertson- Walker 
metric, with the extrinsic curvature and the Ricci tensor to be 



K, 



ij 
R, 



H(t)gij, 



K = 3H(t), 



R 



(40) 
(41) 



where H(t) = a/ a is the Hubble parameter. 

We take the variation of the action (|35|) with respect to N, and have our first equation 
of constraint. 



V9 



-jL^Ky-XK^ + CR + a 



d 6 x = / y/gpd 6 *.. 



(42) 



Here, p is the energy density of the Lifshitz scalar in the universe, and can be written as 



P 



1 SS m 
y/9 SN 



(43) 



where S m is the action of matter field, which can be a Lifshitz scalar (l7| . gauge fieldjij 
or something else. Because of the projectability of the lapse function N(t), we only have 
a spatial-integral constraint here. This is generic for all the Horava-like models with pro- 
jectable lapse function iV(f)(2a|. But, for a homogeneous and isotropic Friedman universe, 
this constraint equation is valid at every point, and the integral can be removed legally. 
Thus we have the first Friedman's equationjlll. Il7 



H 2 



3a(l -3A) 



-p + — e(A) 



a(X) 



(44) 



Since p is the energy density of matter and radiation, a(X) plays the role of "cosmological 
constant". Here, it is a function of A and evolves when A varies as the energy scale changes. 
This implicant dependence may be treated carefully when we are facing problems like the 
evolution of dark energy or the tilt of the power spectrum. But because the dependence 
of A on the cosmic time is unknown, we will neglect this dependence and suppose that in 
the process we are interested in, the change of a(X) is so little that it will not have any 
significant physical effect, and so is H(X). We see from ( 14"4"|) that if the universe is flat 
and dominated by the cosmological constant, for some A greater than 1/3, we must have 
cr(A > 1/3) < 0, which means that we have a positive cosmological constant A > at 
IR, since from (131)1) . £(A > 1/3) is always positive. These two conditions guaranteed the 
positivity of the cosmological constant and H 2 . If the matter/radiation contribution could 
be ignored safely, the homogenous and isotropic solution is a pure de-Sitter spacetime, with 
an exponentially expanding scale factor a(t) cx exp(Ht). 
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The second equation of constraint is obtained by taking the variation of the action with 
respect to the shift vector N i: 

Vi{K lj - XKg lj ) = 0. (45) 

Because the extrinsic curvature is homogeneous in a Friedman universe, as in (j4"0l . oc 
gij, this equation is trivially satisfied for the background evolution. But it will supply a 
perturbative constraint equation up to first order if the perturbations to the background 
metric are under consideration. 

Finally take the variation of action (155]) with respect to g^, we have the equation of 



motion of dynamical degree of freedom. The explicit expression (see [111, [18|, [26J]) is rather 



lengthy and has little to do with our following discussion so we would like not write it here. 



IV. SPECTRUM OF THE GRAVITATIONAL SCALAR 

From the discussion above, we know that the classical evolution of the scale factor in the 
Hofava era is determined by (144"I) . Especially, when cosmological constant is dominant and 
the universe is flat, the evolution is the exponentially expansion like in a de Sitter phase. 
In this section, we will calculate the perturbation of the gravitational field, and study the 
equations of motion of the scalar modes. 

Taken the ADM formalism and the gauge choice (1231) . with some relations derived through 
the two constraint equations, the perturbed action of the gravitational field up to second 
order can be written as 



S (2) 



3 1 - A 



dtd 3 x{ 3aa 3 (l - 3A) 
~(3/3i + 2f3 2 + 8/? 3 )^ 6 V> - 2a£(A)V>#V } . 



2 V 2 



+ QHtptp + 9H 2 ip 2 



(46) 



Now the Hubble parameter is a constant. For convenience we define a conformal time rj with 
dt = adrj and introduce an auxiliary field \ = gm/>. After taking the variation with respect 
to ip, and changing to the momentum space, we have 



where 



is the speed of sound, and 



x "(v) + k H A L rf + c 2 s k 2 - - x (v) = 0. 



3A 



(47) 



(48) 



2tt : 



2tt 



(3 1- A 
a 1 — 3A 



(3A + 2p 2 + 8/3 3 ) 



(49) 



is the characteristic length which denotes the scale where the trans-Planckian effects becomes 
significant. 

This equation can not be solved analytically. However, many efforts has been done to 
deal with this type of Corley-Jacobson dispersion relation [2^, [23|, 24|. The author of 2^ 
have already studied the trans-Planckian physics appearing naturally in Horava-Lifshitz 
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gravitational waves, and find a scale invariant power spectrum in a specific connecting time. 
Here we follow the method of Martin et.al.|22j to investigate the trans-Planckian effects 
of the gravitational scalar. We split the period under consideration into three regions by 
two different characteristic lengths: the trans-Planckian length L and the sound horizon 
c s /H. First, when rj — > — oo, we can set the initial value of the wave function xivd an d its 
derivative x'ivd such that they initially minimize the energy density. This is to satisfy 



x'iVi) 



WHL\rn\ 




±i\ —HL\rji\. 



(50) 
(51) 



Then, in UV region when the physical wavelength of the mode concerned is much smaller 
then the characteristic length, k 6 term is dominant in (j4"7| . 

X '{j V + k 6 H A L A V \uv = 0. (52) 
In a new variable z = k 3 H 2 L 2 \r]\ 3 /3, the solution to this equation can be expressed as 

Xvwijl) = A 1 ^J\q\J 1 / & (z) + A 2 \/\r)\J-i/6(z), (53) 

where Ai and A 2 can be determined by the continuity of x an d x' at the initial time rji, to 
be 



A 1 
A 2 



txHL 



6 sin(7r/6) 

nHL 
6 sin(7r/6) 



/c 3 |r]i| 



■[J-5/e(zi) T^i/e^i)], 



(54) 
(55) 



Note, that in the UV region ^> 1 and z- x ^> 1, we can expand the Bessel's function into 
its asymptotic form when the argument is large, and the coefficients reads 



A l 



A, 



|(^/) PX] > 

^(±i)exp 



71 71 

=F« I Zi H 

1 12 4 

±1 [ Zi 

1 12 4 



Here for simplicity we have defined 



x 



z H , 

12 4' 



7i 

12 



±i 
±i 



7i 

4' 



— e TiXi 
3 



±iyi 



—e 



(56) 
(57) 

(58) 



In the intermediate region, the wavelength of the /c-mode exceeds the characteristic length 
L but still much less then the sound horizon c s /H, i. e. L ^> A ^> c s H~ l . Then we can 
neglect the UV term and the cosmological damping term in f|4T|) . and get an oscillation 
solution of the perturbation as a plane wave. We can deduce the connecting time 77* between 
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UV and intermediate region 



2 

> 



2tt , N r , . 2tt 1 

T a(^) = L„ \v*\ = Wrrmr - (59) 

Here, the emergence of the subscript * of the characteristic length L* is because that L 
also varies with time through the parameter A, which depends on the energy scale by the 
renormalization flow. When we calculate the critical time r)+, all the parameters should be 
the value at 77* including L. In this region, the solution is the usual plane wave 

Xint(v) = Bie ic ^ + B 2 e- ic ° kr >. (60) 

Then we can determine the coefficients B\ and B 2 by connecting Xmt and its derivative with 
those in the UV region. This reads 

+ Y^*\ 1/2 [J-i/M - ik 2 H 2 L 2 \^\ 2 J^)] , (61) 

+ y^\ 1/2 + ik 2 H 2 L 2 \ Vi( \ 2 J 5/6 (z,)] . (62) 

To go further we notice that 

k 2 H 2 L 2 M 2 = 1, (63) 

and generally <C if -1 , so we also have |7^| ^> 1. The Bessel functions can also be 

expanded in its asymptotic form as before, then we get 

B 1 = T -j— e -^v* e ±i*i[ e -iy* _ e T«/8 e -fa*] ) (64) 



B 2 = T ^— e ^kv, e ±i^\ [e w,_ e T^/Q e i^ (g 5 ) 
V 2k 



where x and y are defined in (joojl . 

As the universe expands the wavelength of /c-mode is stretched and becomes larger and 
larger, and finally exceeds the sound horizon. This critical time 77* is determined by 

— =H' lV * l = —k- (66) 

The meaning of the subscript * of the sound speed is similar as before: we require A in the 
definition of c s be its value at 77*. When 77 ^> 77*, the fluctuation is in an IR region. The 
perturbation will freeze out after it exceeds the horizon, so primordial value of the power 



2 Here some subtlety exists when choosing the exact connecting time whether we should let A = L or 
k 6 H i L i r/ i = c? s k 2 . In [23J], the detailed discussion shows that we should choose to avoid the oscillation 
spectrum. But in our case we will show that the former choice will also sweep the oscillation in the final 
result. 
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spectrum observed today can be traced back to its value at and before 77*. The IR solution 
of equation (H7|) is 

Xm = C V 2 + - (67) 
V 

For convenience we only pick the increasing mode 3 , whose coefficient can be determined by 
connecting Xmt and xir at 77*, which reads 



D 



|=Fi fc(?7* - 77*) + - Zi - ^ J. (68) 



Therefore the power spectrum is 

_* Jv i,iO k TT .n. ^.9 H ^ 



27T 2ITI 2vr 2 11 \2tvc 



\W = 7^H 2 \D\ 2 = —— . (69) 



which is obviously scale- invariant, if we neglect the time variation of the horizon at the 
inflationary stage. If so, the slight difference of the horizon-crossing time for different wave- 
lengths will produce different H's, thus different spectra. This is of course the usual case in 
the model slightly breaking the time-translation invariance of the de-Sitter stage. 

To comprehend the significance of ( 1691) , we note that the spectrum is frozen out after the 
mode exceeds the sound horizon. When the exponentially expansion is over and the universe 
recovers usual GR behavior, all the parameters in our original action (1351) recover the IR 
limit value, and specifically, A — > 1 + . But this will never influence the value of the power 
spectrum which is completely determined by the values of parameters at horizon-crossing. 
On the other hand, if A runs to unity much earlier before the wavelength of the fluctuation 
mode exceeds the sound horizon, i. e. c s (r] < 77*) = 0, 4 then the scalar spectrum is divergent, 
showing the breakdown of the treatment. This is just what we expect: the theory recovers 
the general relativity so early that it looks the same as the usual inflationary model and the 
gravitational scalar is not a physical degree of freedom and will not bring observable power 
spectrum any more. 



V. PRIMORDIAL GRAVITATIONAL WAVES 

From the discussion above we see that the scalar mode of the gravitational scalar is scale 
invariant. This is the generic property under the de Sitter background with time translation 
invariance. We will see the same result for the gravitational tensor modes[29j]. We also start 
with the perturbed metric, only to the tensor parts, 

ds 2 = -dt 2 + a 2 (t) (5ij + hij(t, x)) dx i dx j , (70) 

where hy has been already defined in [TUl and satisfies the transverse-traceless conditions 
h\ = 0, dihij = 0. Substituting this metric into the total action, we obtain the tensor action 
of second order, 



= I dtd 3 xa 3 



(71) 



Connect with both the decaying and increasing modes will bring here an inessential factor of order unity, 
as in [? ] . 

This requires c(A) remains finite in the case when A — » 1. 
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The transverse traceless tensor hij can be Fourier transformed by plane waves with wavenum- 
ber k as 

h tJ (t,k)= f ^W^. (72) 

A=R,L ' 



kp A p\ A [25 



where pf^ is the circular polarization tensor which is defined by %k s e rs ip A - 
Here p R = 1, p L = — 1, and are called the right handed mode and the left handed mode 
respectively. We also impose normalization conditions as p* lA p{ B = 5 AB , where p* tA is the 
complex conjugate of p lA . Substituting the expansion into the action ( 1711 . we obtain 



5 S„ 



A=R,L 



f , d 3 k , fa 




2 


/ dt - — —or < — 







4a 6 4a 2 



I , 4 I 2 



(73) 



Using the variable = and the conformal time 77, and rewrite the action by co-moving 
time rj, with a = —1/Hrj in the de Sitter-like space, we have 



v w ( V ) + k b HHY + c 2 k 2 - - v k ( V ) = 0. 



(74) 



where similarly to the treatment on the scalar perturbation before, we define a characteristic 
length related to ultraviolet gravitational waves, I = (j3i/3 / 'a) 1//4 . We see from this definition 
and (I49p . that the characteristic length of scalar and tensor modes may be different. This 
relies on the relative magnitude of different /3's. 

Some discussions parallel to last section will yield the power spectrum of gravitational 
wave as 

, = ^! = fJLf (75) 



2tt 2 



2ttc 



where the subscript t means the quantities are evaluated at the time of horizon-crossing of 
the gravitational waves, i.e. 

IW = — ■ (76) 



A- 



ct/c 



Obviously the power spectrum of tensor mode is scale invariant as well. In [29], the author 
use the same method to connect the solution step by step and calculate the infrared power 
spectrum, but with only two pieces to join and thus more accurate sub-solution in each 
piece: Hankel function in infrared region. However, after taking the correct connecting time 
V* m (]59j) as we do before, the dependence on the "cutoff energy scale" there also vanishes. 
Now we can calculate the tensor-to-scalar ratio 5 , 



ct 



1 - A* e(AQ 

l-3A*£(A t )' 



(77) 



5 Note that the tensor-to-scalar ratio is usually defined to be ^h/^Stp, by the power spectrum of the 
perturbations of inflaton field in ordinary inflationary models. And in super-horizon scales, ^s<)> is of the 
order as 3?h divided by slow-roll parameter e. Because we have not placed any scalar field here, we define 
the tensor-to-scalar ratio as S?h divided by just the spectrum of gravitational scalar 
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This shows that the ratio is only determined by the speed of sound/light at the sound- 
/hubble-horizon-crossing. All the dependence on the characteristic length of the ultraviolet 
behavior do not appear in the ratio. We may try to estimate this ratio by the assumption 
that £ varies slowly to 1 when the k-mode we are interested in crosses the horizons. Then 
we can neglect the £ term in (1771) and have only the dependence of A at sound horizon 
crossing. Since r«l from the observations, A* = 1 + 2r must be very close to 1 + , which is 
in consistency with our former assumption. 



VI. CONCLUSION AND DISCUSSION 

In this paper, we clarified several issues in the Hofava-Lifshitz gravity. We first showed 
that the strong coupling issue may not be so serious as argued in the literature before. The 
basic point is that the diffeomorphism is only a good approximation even at IR. Taking 
into account of the projectability condition, the usual Stuckelberg trick could not be applied 
naively. From our discussion, it seems that the extra dynamical scalar degree of freedom 
could be decoupled naturally. 

However, the theory may suffer from other pathologies. One concern is on the existence 
of the ghost excitation. We showed that as the perturbations around the flat spacetime, the 
scalar perturbation around the flat FRW universe could be a ghost in the parameter region 
| < A < 1. The presence of the ghost mode is a serious challenge to the theory. We tried 
to avoid the dangerous parameter region by mildly modifying the Hofava-Lifshitz action. 
We kept only the most UV sensitive and IR sensitive terms. We discussed the classical 
evolution and the power spectra of scalar and tensor perturbations. We obtained scale 
invariant spectrum if the Hubble constant H does not change. We also calculated the 
tensor-scalar ratio, and found it could be small under reasonable condition. 

The nature of the power spectra studied in this paper is purely gravitational. In par- 
ticular, in the language of orthodox cosmology, the scalar perturbation is expected to set 
up the initial conditions and seed the anisotropy of large scalar structure in our universe. 
Some work has been done to reveal the evolution of perturbations after inflation in Hofava- 



Lifshitz gravity 32] . After inflation ends, this gravitational perturbation must be converted 
into CMB anisotropy and matter inhomogeneity through some post-inflation evolotions. But 
still we do not know yet how to couple the gravitational scalar mode with, for instance, the 
radiation. This is an interesting issue, which we would like to study in future. Recently, 
an interesting paper on the Hofava-Lifshitz universe with single scalar field discussed the 
curvature perturbation CjHj]. 

One essential issue in Hofava-Lifshitz gravity is on its RG flow. In [27J , it has been shown 
that in Lifshitz-like scalar field theory, the RG flow may not lead the theory to the fixed point 
we want. Considering the numbers of the parameters in modified Hofava-Lifshitz gravity, 
this raise the concern if the theory can flow to IR fixed point A = 1. Moreover, the details of 
RG flow can tell us if we can avoid the dangerous region, where the ghost excitation appears, 
even we start from a safe region. Furthermore, RG flow may closely related to the physics 
in the inflationary era. It is not clear whether RG flow of the theory runs to its IR fixed 
point before the inflationary era. If it did, then the gravitational scalar is not dynamical 
and has nothing to do with inflation. Even if the energy scale to reach IR limit is lower than 
the inflation era, there is an important question to answer: did A vary significantly in the 
inflationary era? The variation of A may tilt the power spectra and has interesting physical 
implications. In any case, the behavior of the Hofava-Lifshitz gravity theory under RG flow 



15 



deserves careful investigations. 



Acknowledgments 

The work was partially supported by NSFC Grant No. 10535060, 10775002, 10975005 
and RFDP. BC would like to thank R.G. Cai, Miao Li and Jian-xin Lu for stimulating 
discussions. BC also thank ICTS-USTC for hospitality, where the project was inspired. SP 
is grateful to Yi-Fu Cai and Tower Wang for useful directions. 



[1] P. Horava, "Membranes at Quantum Criticality," JHEP 0903, 020 (2009) |arXiv:0812.4287l 
[hep-th]]. 

[2] P. Horava, "Quantum Gravity at a Lifshitz Point," Phys. Rev. D 79, 084008 (2009) 

|arXiv:0901.3775l [hep-th]]. 
[3] E.M. Lifshitz, "On the Theory of Second-Order Phase Transitions I &; II", Zh. Eksp. Teor. 

Fiz 11 (1941)255 & 269. 
[4] B. Chen and Q. G. Huang, "Field Theory at a Lifshitz Point," !arXiv:0904.4565l [hep-th]. 
[5] R. G. Cai, B. Hu and H. B. Zhang, "Dynamical Scalar Degree of Freedom in Horava-Lifshitz 

Gravity," larXiv:0905.0255l [hep-th]. 
[6] B. Chen, S. Pi and J. Z. Tang, "Scale Invariant Power Spectrum in Hofava-Lifshitz Cosmology 

without Matter," larXiv:0905.2300l [hep-th]. 
[7] C. Charmousis, G. Niz, A. Padilla and P. M. Saffin, "Strong coupling in Horava gravity," 

larXiv:0905.2579l [hep-th]. 

D. Bias, O. Pujolas and S. Sibiryakov, "On the Extra Mode and Inconsistency of Horava 

Gravity," larXiv: 0906 .30461 [hep-th]. 
[8] M. Fierz and W. Pauli, "On relativistic wave equations for particles of arbitrary spin in an 

electromagnetic field," Proc. Roy. Soc. Lond. A 173, 211 (1939). 
[9] N. Arkani-Hamed, H. Georgi and M. D. Schwartz, "Effective field theory for massive gravitons 

and gravity in theory space," Annals Phys. 305, 96 (2003) | arXiv:hep-th/0210184| . 
[10] H. van Dam and M. J. G. Veltman, "Massive And Massless Yang-Mills And Gravitational 

Fields," Nucl. Phys. B 22, 397 (1970). 
[11] H. Lu, J. Mei and C. N. Pope, "Solutions to Horava Gravity," larXiv:0904.T59~5l [hep-th]. 
[12] A. Kehagias and K. Sfetsos, "The black hole and FRW geometries of non-relativistic gravity," 

Phys. Lett. B 678, 123 (2009) [arXiv: 0905 .04771 [hep-th]]. 
[13] J. Z. Tang and B. Chen, "Static Spherically Symmetric Solutions to modified Horava-Lifshitz 

Gravity with Projectability Condition," larXiv:09093l~27l [hep-th]. 
[14] T. Harko, Z. Kovacs and F. S. N. Lobo, "Testing Hofava-Lifshitz gravity using thin accretion 

disk properties," Phys. Rev. D 80, 044021 (2009) |arXiv:0907. 14491 [gr-qc]]. 

T. Harko, Z. Kovacs and F. S. N. Lobo, "Solar system tests of Hofava-Lifshitz gravity," 

larXiv: 0908 .2874 [gr-qc]. 

L. Iorio and M. L. Ruggiero, "Horava-Lifshitz gravity and Solar System orbital motions," 
larXiv:0909.2562l [gr-qc] . 
[15] M. Li and Y. Pang, "A Trouble with Hofava-Lifshitz Gravity," larXiv:0905."2~75"T1 [hep-th]. 



16 



[16] Y. W. Kim, H. W. Lee and Y. S. Myung, "Nonpropagation of scalar in the deformed Hofava- 
Lifshitz gravity," larXiv:0905.3423l [hep-th]. 

[17] G. Calcagni, "Cosmology of the Lifshitz universe," larXiv:0904.0829l [hep-th]. 

[18] E. Kiritsis and G. Kofinas, "Horava-Lifshitz Cosmology." llarXiv:0904.1334l [hep-th]. 

[19] S. Mukohyama, "Scale-invariant cosmological perturbations from Horava-Lifshitz gravity with- 
out inflation," larXiv:0904.2l90l [hep-th]. 

[20] R. Brandenberger, "Matter Bounce in Horava-Lifshitz Cosmology," arXiv:0904.2835 [hep-th]. 
Y. S. Piao, "Primordial Perturbation in Horava-Lifshitz Cosmology," arXiv:0904.4117 [hep- 
th]. X. Gao, "Cosmological Perturbations and Non-Gaussianities in Horava-Lifshitz Gravity," 
arXiv:0904.4187 [hep-th]. T. Sotiriou, M. Visser, S. Weinfurtner, "Phenomenologically viable 
Lorentz- violating quantum gravity", arXiv:0904.4464 [hep-th]. S. Mukohyama, K. Nakayama, 
F. Takahashi and S. Yokoyama, "Phenomenological Aspects of Horava-Lifshitz Cosmology", 
arXiv:0905.0055 [hep-th]. S.Kalyana Rama, "Anisotropic Cosmology and (Super) Stiff Matter 
in Hofava's Gravity Theory", arXiv:0905.0700 [hep-th]. Alex Kehagias, Konstadinos Sfetsos, 
" The black hole and FRW geometries of non-relativistic gravity", arXiv:0905.0477 [hep-th]. 
E. N. Saridakis, "Horava-Lifshitz Dark Energy," arXiv:0905.3"532l [hep-th]. X. Gao, Y. Wang, 
R. Brandenberger and A. Riotto, "Cosmological Perturbations in Horava-Lifshitz Gravity," 
arXiv:0905.3821 [hep-th]. M. Minamitsuji, "Classification of cosmology with arbitrary mat- 
ter in the Horava-Lifshitz theory," larXiv:0 905,3892 [astro-ph.CO]. A. Wang and Y. Wu, 
"Thermodynamics and classification of cosmological models in the Horava-Lifshitz theory 
of gravity," arXiv:0905.4117 [hep-th]. S. Nojiri and S. D. Odintsov, "Covariant Horava-like 
renormalizable gravity and its FRW cosmology," arXiv:0905.4213 [hep-th]. M. i. Park, "The 
Black Hole and Cosmological Solutions in IR modified Horava Gravity," arXiv:0905.4480 [hep- 
th]. Y. F. Cai and X. Zhang, "Primordial perturbation with a modified dispersion relation," 
larXiv: 0906 .334T1 [astro-ph.CO]. M. i. Park, "A Test of Horava Gravity: The Dark Energy," 
arXiv:0906.4275 [hep-th]. Y. Lu and Y. S. Piao, "Scale Invariance from Modified Dispersion 
Relations," arXiv:0907.3982 [hep-th]. C. Appignani, R. Casadio and S. Shankaranarayanan, 
"The Cosmological Constant and Horava-Lifshitz Gravity," arXiv:0907.3121 [hep-th]. 

[21] C. Bogdanos and E. N. Saridakis, "Perturbative instabilities in Horava gravity," 
larXiv: 0907. 16361 [hep-th]. 

[22] J. Martin and R. H. Brandenberger, "The Trans-Planckian Problem of Inflationary Cosmol- 
ogy", Phys. Rev. D 63, 123501(2001) [arXiv:hep-th/0005209] . 

[23] J. Martin and R. H. Brandenberger, "The Corley-Jacobson dispersion relation and trans- 

Planckian inflation," Phys. Rev. D 65, 103514 (2002) |arXiv:hep-th/0201189| . 
[24] R. H. Brandenberger and J. Martin, Mod. Phys. Lett. A 16, 999 (2001) 

[arXiv:astro-ph/0005432] . 
[25] T. Takahashi and J. Soda, "Chiral Primordial Gravitational Waves from a Lifshitz Point," 

larXiv:0904.0554l [hep-th]. 
[26] T. P. Sotiriou, M. Visser and S. Weinfurtner, "Quantum gravity without Lorentz invariance," 

larXiv:0905.2798l [hep-th]. 
[27] R. Iengo, J. G. Russo and M. Serone, "Renormalization group in Lifshitz-type theories," 

larXiv:0906.3477l [hep-th]. 
[28] G. Calcagni, "Detailed balance in Horava-Lifshitz gravity," arXiv:0905.3740 [hep-th]. 
[29] S. Koh, "Relic gravitational wave spectrum, the trans-Planckian physics and Horava-Lifshitz 

gravity," larXiv:0907.0850l [hep-th]. 



17 



[30] K. Yamamoto, T. Kobayashi and G. Nakamura, "Breaking the scale invariance of the pri- 
mordial power spectrum in Horava-Lifshitz Cosmology," Phys. Rev. D 80, 063514 (2009) 
|arXiv:0907.1549 [astro-ph.CO]]. 

[31] A. Wang and R. Maartens, "Linear perturbations of cosmological models in the Horava-Lifshitz 
theory of gravity without detailed balance," arXiv:0907.1748 [hep-th]. 

[32] T. Kobayashi, Y. Urakawa and M. Yamaguchi, "Large scale evolution of the curvature per- 
turbation in Horava-Lifshitz cosmology," arXiv:0908.1005 [astro-ph.CO]. 

[33] S. Carloni, E. Elizalde and P. J. Silva, "An analysis of the phase space of Horava-Lifshitz 
cosmologies," larXiv:0909.221~9l [hep-th]. 

[34] A. Wang, D. Wands and R. Maartens, "Scalar field perturbations in Horava-Lifshitz cosmol- 
ogy," larXiv:0909.5T67l [hep-th]. 



18 



